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The modified _F(_R)-scalar-Gauss-Bonnet gravity is proposed as dark energy model. The recon- 
struction program for such theory is developed. It is explicitly demonstrated that the known classical 
' universe expansion history (deceleration epoch, transition to acceleration and effective quintessence, 

I phantom or cosmological constant era) may naturally occur in such unified theory for some (re- 

, constructed) classes of scalar potentials. Gauss-Bonnet assisted dark energy is also proposed. The 

^ ' possibility of cosmic acceleration is studied there. 
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I. INTRODUCTION 

^ i' It is quite possible that dark energy is some manifestation of the unknown gravitational physics. The search 
D \ of the realistic gravitational alternative for dark energy (for a recent review, see continues. According to this 
approach, usual General Relativity which was valid at deceleration epoch should be modified by some gravitational 
[ terms which became relevant at current, accelerating universe when curvature is getting smaller. Among the most 
■ popular modified gravities which may successfully describe the cosmic speed-up is F{R) gravity. Very simple versions 
> ' of such theory like 1/i? and 1/R + H may lead to the effective quintessence/phantom late-time universe. 
1 Moreover, positive/negative powers of curvature in the effective gravitational action may have stringy/M-theory 
' origin [J. Recently, big number of works was devoted to the study of late-time cosmology and solar system tests 
' checks in modified F{R) gravity 0, Q. While it is not easy to satisfy all known solar system tests at once, it is 
possible to construct the explicit models 0,9 which describe the realistic universe expansion history (radiation/matter 
] dominance, transition to acceleration and accelerating era). 

Another theory proposed as gravitational dark energy is scalar-Gauss-Bonnet gravity [9l] which is closely related 
with low-energy string effective action. The late-time universe evolution and comparison with astrophysical data in 
such model was discussed in refs. [ol. [lOl [Til fl^ . The possibility to extend such consideration to third order (curvature 
^ cubic) terms in low-e nerg y string effective action exists [l3| . Moreover, one can develop the reconstruction program 
k> for such theories (see |l4l |. for a review). It has been demostrated that some scalar-Gauss-Bonnet gravities may 
; I ■ be compatible with the known classical history of the universe expansion. 

\ In the present paper, we propose the unified i^(i?)-scalar-Gauss-Bonnet gravity as dark energy model. The recon- 
struction program for such model is explicitly developed. It is shown that it is easier to realize the known, classical 
universe history (deceleration, transition to acceleration and cosmic acceleration with effective w close to —1) in such 
a model than in F{R) or scalar-Gauss-Bonnet gravity. Several different late-time regimes (effective quintessence, 
phantom or cosmological constant) are investigated. Moreover, it is indicated that it is possible to pass Solar System 
tests in such a unified model. It is shown also that Gauss-Bonnet term may play an important role for other class of 
gravitational models where it couples with scalar kinetic term. The possibility of cosmic acceleration in such model 
is also demonstrated. 
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II. COSMIC ACCELERATION IN F(i?)-SCALAR-GAUSS-BONNET GRAVITY 

In this section, let us introduce F(i?)-scalar-Gauss-Bonnet gravity whose action is given by 



_R , r\ 

2k, 



(1) 



Here F(R) is a function of the scalar curvature R, Lmattor is matter Lagrangian density, and G is GB invariant: 
G = i?2 - 4Rf,,Rf''' + R^.p^Rf^'^P''. 

In ([1]), with F(R) — rj = 1 corresponds to string-inspired scalar-Gauss-Bonnet gravity which has been proposed 
as dark energy model [9|| (it may be applied also for resolution of the initial singularity in early universe [111) and 
77 = to /(G) gravity [17[- Note that using trick of ref. 3] one can delete F{R) term, at the price of the introduction 
of second scalar. However, in such formulation the coupling (including derivatives) of new scalar with GB sector 
appears. 

In the following, the metric is assumed to be in the spatially-flat FRW form: 



The FRW field equations look as 



ds' = -dt^ + a{tf {dx'f . (2) 



= - F{R) + 6 + (^-^ + F'iR)^ - 36 (aH^H + HH^ F" 



(R) 



^ ^+FiR)-2{H + 3H^) 

-f48 (iH^H + + 2HH^ F"{R) + 72 (iHH + F"'{R) + \<f^ ~ ^(0) 

dt^ dt at 

and the scalar field equation is 

O = r/(0 + 3i7</)) +T/"'(</.)+e'(0)G . (5) 

Here R = 12H^ + QH and G = 24 [hH^ + iJ*) 

We now consider the perfect fluids with constant equation of state (EoS) parameters Wi = pi/ pi as the matter. The 
energy density is pi — pi^a'^'^^'^'^^'^ with a constant p^Q. Let us parametrize the model with two proper functions /(</>) 
and g{t) as follows (compare with [l^ for F = case), 

+36 (4.9' {f{^)f g" {fm + g' (/(0)) g'" (/(0))) F" {R) - 3g' (/(0)) e^^f^^^^U{ct>) 



1 f^A /'('/'i)e'^(^^*^»t^(0i) 
#1 9 

9' (/(0i))' 



C/(0) ^ y'*d0i/'(0i)e-3(/(^^)) (^4g"(/(0i)) (^^ + F'(i?)^ 

+ 12 (4.9' {f{cb,)f g" (/(0i)) + 45" (/(0i))2 + 5g' (/(0i)) g'" (/(0i)) ) F"{R) 
+72 (4g' (/(<^i)) 5" (/(0i)) + 5"' (/(0i))) ^^'"(i?) 



?7 



/'( 



^2 



+ 5^(1 + z«,)P«ao '''-*-"")e-3(i+-)s(/(^^)) 



R ^ l2g'{f{cj>)f + Qg"{f{<P)) . (6) 
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Here, the equations have the following solution: 

a = aoe^(*) (H^g'it)) , <j> = r\t) . (7) 

This may be considered as reconstruction of above modified gravity from known universe history expansion. 

In case ry = in ([T]) , we can freely redefine the scalar field as —>(/>' = </>'(</>). If depends on the time t as = 
one can redefine the scalar field properly and identify the scalar field as t, (p = t, that is /(0) = </>. Then Eq.® can 
be simplified as (for the case F — compare with [14]) 

+36 (45' {cf>f g" (0) + g' (0) g'" (0)) F" {R) - 3g' {^) ef(^'t/(<^) , 



9' 

U{4>) ^ J dcb^e-^M (^Ag"{cf>,)(^^ + F'{R) 

+ 12 (45' g" (0i) + Ag" {cj>,f + 5.g' (0i) g'" (0i)) F"{R) 

+72(45'(</.i)g"(</.i)+5"' {<t>{)) F'" {R) 

+ ^(l + u;,)p,o«o'''+""'e-3(i+-).(0O^ , 

R = I2g'{4>f+&g"{<f) . (8) 
As some example oi F{R), we consider 0] 

J^(^) = • (9) 

The original model has two solutions corresponding to deSitter universe, where R is constant, and asymptotic universe 
with effective w = —2/3. In the model ([S]), one can realize any time developement of the universe. For example, by 
choosing ^ and 



and /(0) to be properly defined, we obtain 



g{t) = Hot + Hi\n[^] , (10) 



H{t) =Ho + ^ . (11) 



When t is small, H pT|) behaves as that in universe with perfect fluid with w = —1 + 2/3Hi and when t is large, H 
behaves as in deSitter space, where _ff is a constant. Then if we choose Hi = 2/3, we find that before the acceleration 
epoch, the universe behaves as matter dominated one with w = 0. After that, it enters to acceleration phase. In the 
original model 0], it was difficult to realize the matter dominated phase. It is easy to see that matter dominance phase 
(with subsequent acceleration) can be easily realized by adding the Gauss-Bonnet term (see Appendix for explicit 
form of scalar potentials). 

As another example, we consider a model with ^ but 

g{t)=Holnl^Hiln(*-^) , (12) 



to \ to 

which gives 

Here Ho, Hi, and to are positive constants. When t is small, H (|13p behaves in a way corresponding to the perfect 
fluid with w = —I + 2/SHo. Then if we choose Ho = 2/3, we find the matter dominated universe. On the other 
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hand, when t ^ to is large, H behaves as that in the phantom universe with w = — 1 — 2/3-ffi < — 1 and there will 
appear a big rip singularity at t — tQ. The three- year WMAP data are analyzed in Ref. which show that the 
combined analysis of WMAP with the supernova Legacy survey (SNLS) constrains the dark energy equation of state 
wd e pushing it clearly towards the cosmological constant value. The marginalized best fit values of the equation of 
state parameter at 68% confidence level are given by —1.14 < wde < —0.93, which corresponds to Hi > 10.7 as Hi is 
positive. In case one takes as a prior that the universe is flat, the combined data gives —1.06 < wde < —0.90, which 
corresponds to Hi > 25.0. Therefore the possibility that wde < —1 has not been excluded. As clear from ^ or ([5]), 
the expressions of V{(j)) and ^i(0) depend on the explicit form of F{R), say /i in ([5]), and the time-development of 
the universe g{t), now Hq and Hi, (and /(</>) in case of Then the form of F{R) is irrelevant to the WMAP data 
but only g{t) is relevant as long as we use the expressions in ([5]) or ([5]). 
One may also consider a model 

ff(0=i?oln^ + (i?i-i/o)ln(^^) , (14) 
with constants Hq, Hi, and Iq. Then we obtain 

When t is small, H again behaves as _ff ^ Hg/t corresponding to the universe filled with perfect fluid w = —l + 2/3Ho- 
On the other hand, when t ~ to is large, H behaves ds H ^ Hi/t corresponding to the fluid with w ^ — 1 + 2/3Hi. 
Therefore with the choice Hq = 2/3, we find the matter dominated universe in the early universe and with the choice 
Hi > 1, we obtain acceleratedly expanding universe, where w < —1/3. The constraint from only the WMAP data 
indicates Hi > 21.4 and that from the combined data indicates Hi > 15.0. 
One more example is ACDM-type cosmology: 



3(1 + w) V \ 21 ' "'J J 3 



Here / is the length scale given by cosmological constant I ~ (lO^^'^ eV) ^ and Iq is a constant. The time- development 
of the universe given by g{t) pB]) can be realized in the usual Einstein gravity with a cosmological constant A and 
cold dark matter (CDM), which could be regarded as dust. Then in the present formulation, by using V{4>) and £,i{(j)) 
in ([6]) or ([8]), ACDM-type cosmology can be realized without introducing CDM as a matter. 
One may also consider more general choice of F{R) as in [H 

F{R)^-^+/3R\ (17) 



or 



F{R)^^^+PR^" . (18) 



Even in this case, we can realize any time development of the universe. For example, if we choose g{t) as in (|10p . 
H pT|) behaves as that in universe with perfect fluid with w — —1 + 2/3Hi when t is small (especially by choosing 
Hi = 2/3, the matter dominated phase occurs) and H behaves as in deSitter space when t is large. On the other 
hand, if we choose g{t) as in ()12|) . we obtain a model showing the transition from the matter dominated phase to the 
phantom phase. Moreover, with the choice of (I15p . the transition from the matter dominated phase to the quintessence 
phase could be realized. 

Let us now consider the string- inspired modelj^, where 



2<p 2(/» 



V = Voe-^ , = ^oe^ , (19) 

with constant parameters Vq, ^q, and <j>Q (with F{R) given by pT]) ). Different from the model ([6]), it is not straight- 
forward to solve the FRW equations (jl]) and the scalar field equation in this model and to find the behavior of 
the universe. In case |-^| , |/3i?^| ^ |^|, the cosmology [§| could be reproduced, that is, if we choose. 



6 

tl «2(l-H/lo) 



mohl _ ^ fuo- ) . (20) 
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we obtain a solution 

H^^, = 0olnf, (21) 

when /zo > or 

H = = (^oln^^, (22) 

when ho < and a constant ti. We should also note that there could also be deSitter solution: H — Hq , cj) — pq. 
Then the FRW equation ([3]) and the scalar field equation (O have the following simple forms: 

= -4^0 + + l^oe-'P"/^" , = Voe-^P°/'f"> - ^oe'^"/^" . (23) 

K Oil Q 

When Vo and are positive, by combining two equations in (j23p . we obtain 

= -(A-2V6T^)i/„^ + ^. (24) 
If a is positive and -4- > 2-76Vb£o, there can be a solution iJn = -j-^t — " , ^ \ , which gives e^P°/'f'° = 

8^-^75-— 2\/6VnCoj 

^ — 2A/6Vb^o) 3^;^- Therefore, in general, there could occur deSitter solution as late-time universe. 

Hence, we demonstrated that in unified _F(i?)-scalar-GB gravity which may be considered as string-inspired theory 
one may have the number of dark energy scenarios of different type (effective quintessence, effective cosmological 
constant or effective phantom). The corresponding classes of scalar potentials may be easily constructed (see an 
explicit example of such construction in the Appendix). It is interesting that in such unified model it is easier 
to realize the known classical universe history ( radiation/matter dominance, decceleration-acceleration transition, 
cosmic acceleration) than in pure F{R) gravity 0,13] or in pure scalar-GB gravity [3, llBl- There is no problem to 
take into account usual matter in such consideration. In this case it is even easier to reconstruct the decceleration 
(radiation/matter dominance) phase before the acceleration epoch. However, the explicit form of corresponding 
potentials is quite complicated. 



III. GAUSS-BONNET GRAVITY ASSISTED DARK ENERGY 



In this section we study modified gravity model motivated by gravity assisted dark energy |19l | where scalar kinetic 
term couples with the function of Gauss-Bonnet invariant. It is different from the model of previous section, but GB 
term plays again the important role here. The starting action is: 



S 



,R-f{G)Ld 



(25) 



where Ld = ~^g^^'^ di^i4)d^(t>. Taking the same FRW metric and assuming </) only depends on time coordinate t, (f) = (f>{t), 
we find the solution of scalar field equation: 



<A = /(G) 



-1^-3 



Q ■ 



Here g is a constant. 

The FRW field equation is found to be: 



q 



i + 12i/2(i/ + 7ii-2)/'(G) 



/(G) 



2H^) - 12H 



/(G) 



(26) 



0. (27) 



Here /' = df{G)/dG. The following choice of /(G) is convenient: /(G) — • Here a (non-integer) and ^'^ are 

constants. Then in order that /(G) is real, we only consider the case that G has a definite sign and when G is positive 
(negative), fi'^ > (/i^ < 0). One can see that for a = — 1, FRW equation takes the very simple form: 



k2 



288g2ij4 



(28) 
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When fi"^ > 0, the equation has only trivial solution H = 0. When fi* < 0, the following solution exists: 



1/3 



(29) 



which gives H = l/3t and G = -16/27t''. 

Now for general a, by using a constant x, we assume 



a = a^t 



(30) 



Then it follows 



_ 2a+l 
3 ' 



\(2a+l)(a-l) / 



(ct''-3a-l) 



Since G = i^i^fH^Kfi; 



i)r 



— , we find /i^ > when 
the universe accelerates 



(a-l) 

a > 1 or a < -1/2 and < when -1/2 < a < 1. Since a = |ao(2Q; + l)(a 
(a > 0) if a > 1 or a < —1/2. We should also note that 81/x'*/(2q! + l)(a — 1) is positive by assumption. Since Oq > 0, 
one gets (3 — v'l3)/2 < a < 1 or a > (3 + v'l3)/2. Hence, if (3 — \/T3)/2 < a < 1, the expansion of the universe is 
decelerating but if a > (3 + vT3)/2, the expansion is accelerating. 

The effective equation of state parameter in our case \s w = . If a < — i the effective phantom era occurs, 
while if a > 1 the effective quintessence era emerges. 

Let us add the perfect fluid p = (7 — l)p, where 7 g (0, 2), to our system. The FRW equation may be presented in 
the following form: 



ij2 

— = P + PG 



where 
PG = 



a{a + 1)HH^ + 2{2a^ + 6a + 1)HH^ + 2(a + 1) ( a + ^ ) H^H^ + (1 + 7a)H*^ 



(31) 



(32) 



The energy density p satisfies the equation p = —3'yHp. From (|3ip . it follows p — H^/k^ — pQ. Substituting it into 
the conservation law, one gets: 2-^pi — Pg + •^1^ — S^Hpc ~ 0. By assuming we find a — aot^^— as in the 
case without the perfect fluid and p = po^~^^^^^"^- If p oc pc, we find 7(1 + 2a) = 2. For the case that matter is 
dust with 7 = 1, if a > ^, the matter could dominate in the early universe {t — s- 0) and pc dominates in the late 
universe t +00. Hence, if a > (3 + \/T3)/2, there could occur the transition from matter dominated phase to the 
accelerated expansion. It is interesting to note that de Sitter solution is impossible in this scenario with only pc. 

Thus, we demonstrated the possibility to realize the Gauss-Bonnet gravity assisted dark energy with late-time cosmic 
acceleration. The decceleration era may naturally emerge in such scenario before the acceleration with subsequent 
transition to acceleration. 



IV. DISCUSSION 



In summary, we studied late-time, dark energy era in i^(_R)-scalar-Gauss-Bonnet gravity. The reconstruction method 
for such model is developed. It is shown that it may be reconstructed from the known universe expansion history so 
that for some class of scalar potentials the radiation/matter dominance is realized subsequently transiting to cosmic 
acceleration. Moreover, any type of cosmic speed-up: the effective quintessence, phantom or A-CDM era may occur 
after deceleration epoch. It is remarkable that it is easier to achieve the deceleration era transiting to dark era than 
in F{R) or scalar-Gauss-Bonnet gravity. 

Gauss-Bonnet gravity assisted dark energy is proposed. It is shown that cosmic acceleration may naturally occur in 
such theory as well. It is interesting that adding the scalar potential to kinetic term in such theory one can construct 
the Gauss-Bonnet induced model for the dynamical origin of the effective cosmological constant in close analogy with 
ref.0. 

Let us discuss now solar system tests for our mod_el. It has been argued recently that F(i?)-gravity does not pass 
the solar system tests [23]. Recently, however, in j2l|, the conditions that even pure i^(i?)-gravity (of special form) 
could satisfy the solar system and cosmological tests are derived. The conditions are: 

lim F{R) = const , lim F{R) = . (33) 
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An explicit example of such function (with specific values of parameters) is given by 

Hence, we can choose F{R) as p4p and then apply the formulations ([6|) (for 77 = 1 case) and ([8]) (for 77 = case) 
and therefore we can realize arbitrary (accelerating) cosmology by properly choosing V{(j)) and f (</>). In this case, 
our theory passes solar system tests. For F(G')-gravity, which corresponds to 77 = case, it has been shown that the 
scalar field (j) does not propagate and the only propagating mode is graviton, at least in the de Sitter background. The 
model proposed in this paper is the hybrid of i^(_R)-gravity and the scalar-Gauss-Bonnet or _F(G')-gravity. Hence, the 
above arguments indicate that cosmological and solar system tests could be satisfied at once if we choose F{R) (as 
above), V{(j)), and C(0): properly. Of course, to satisfy all cosmological tests some tuning of parameters of our theory 
may be necessary what should be checked explicitly for any particular model. 
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APPENDIX A 



In this Appendix for the simple cosmology case ^ and ([TO]) , we present the explicit form of V{4)) and ^1. 
For the F(i?)-scalar-Gauss-Bonnct model (77 = 1 in ([9])), we further choose /(0) — fo4>- Then V{4)) and ([6]) have 
the following form: 



2?^- 



12(i^o+#i:) - 



- 6 I ( i^o + ^ I - 



+72 -4 ffo 



6gi 



+ 2 Lff + 



Hi 

Hi\ Hi\ 



1 

2^ 



l2(Ho + ^\ - 



6H1 



fo<t> J JfF j J 



12 Ho + ^] -W^ 



fo4>J fS'l^^ 



Hi 



ui<p) 



6(0) - 



C/(0) 



d<Pi 



iM 



I 



-Hi 



-4 



Hi 



V 



1 

2^ 



12 H, + ^ 



6H1 



/o0i J f'sVi J J 



-24 \ -A\H, 



HiY Hi 



-F432 l-ilH 



MJ fir fl<t>\ 
Hi\ Hi Hi 



10 H 



H, \ H, 



20 



MiJ fi<l>' 



_ 6H1 



MJ fir fir 



12(i/o + #i) - 



6gi 
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/ j:- , \ -3{l+Wi)Hi\ 

+ jl + E(l + ^.)p.oao ''^+-)e-3(i+-)^°/"* [^l^j j . (Al) 

This defines the class of scalar potentials which lead to above cosmological solution. For /(G) gravity with f] ^ 0, the 
similar form of V{(f)) and f (0) may be derived. 



[1] S. Nojiri and S. D. Odintsov, Int. J. Geom. Meth. Mod. Phys. 4, 115 (2007) arXiv:hep-th/0601213 . 

[2] S. Capozziello, Int. J. Mod. Phys. D 11, 483 (2002); S. Capozziello, S. Carloni and A. Troisi, arXiv: astro-ph /030304T] 

S. M. Carroll, V. Duvvuri, M. Trodden and S. Turner, Phys. Rev. D 70 (2004) 043528. 
[3] S. Nojiri and S. D. Odintsov, Phys. Rev. D 68, 123512 (2003) arXiv:hep-th/0307288^. 
[4] S. Nojiri and S. D. Odin tsov, Phys. Lett. B 576, 5 (2003) arXiv:hep-th/030707Tp 

[5] F. Faraoni, |arXiv:g r-qc/0607116 arXiv:gr-qc/0511094 ; M. Ruggiero and L. lorio, 'arXiv:gr-qc/0607093 ; A. Cruz- 
Dombriz and A. Dobado, arXiv:gr-qc/0607118 ; N. Poplawski, arXiv:gr-qc/0607124; A. Brookfield, C. van de Bruck 
and L. Hall, arXiv:hep-th/0608015 , Y. Song, W. Hu and I. Sawicki, arXiv:astro-ph/0610532 ; B. Li, K. Chan and 
M. Chu, arXiv;astro-ph/0610794, X. Jin, D. Liu and X^Li, arXiv:astro-ph/0610854 ; T. Sotiriou and S. Liberati, 

[arXiv:gr-qc/0604006; T. Sotiriou, arXiv:gr-qc/0 611107| |arXiv: gr-qc/0604028 R. Bean, D. Bernat, L. Pogosian, A. Sil- 
vestri and M. Trodden, arXiv:astro-ph/0611321 ; I. Navarro and K. Van Acoleyen, arXiv:gr-qc/0611127; A. Bustelo 
and D. Barraco, arXiv:gr-qc/0611149 ; G. Olmo, arX iv:gr-qc/0612047 | T. Saidov and A. Zhuk, arXiv:hep-th/0612227| 

|arXiv:hep-th/0604131 ; F. Briscese, E. Elizalde, S. Nojiri and S7 D. Odintsov, Phys. Lett. B 646, 105 (2007) 
larXiv:hep-th/0612220 ; L. Amendola, R. Gannouji, D. Polarski and S. Tsujikawa, arXiv:gr-qc/0612180 ; S. Baghram, 
M. Farhang and S. Rahvar, arXiv;astro-ph/0701013 , D. Bazeia, B. Carneiro da Cunha, R. Menezes and A. Petrov, 

,arXiv :hep-th/0701106 P. Zhang, arXiv:astro-ph/0701662 ; B. Li and J. Barro w, arXiv:gr-qc /0701111 S. Blud- 
man, arXiv:astro-ph/0702085 , I. Sa wicki and W. Hu, |arXiv:astro-ph/0702278| T. Rado r | ar Xiv:hep-th/0702081 

|arXiv:hep-th/0 701267t L. Sokolowski, 'arXiv:gr-qc/0702097( V. Faraoni, 'arXiv:gr-qc/0703044' 'arX iv:astro-ph /0610734' 
S? Rahvar and Y. Sobouti, arXiv:astro-ph/0704.0680; O. Bertolami, C. Boehmer, T. Harko and F. Lobo, arXiv:gr- 
qc/0704.1733 . 

[6] S. Nojiri and S. Odintsov, Gen. Rel. Grav. 36, 1765 (2004) (arXiv:hep-th/0308176' ; P. Wang and X. Meng, 
|arXiv:astro-ph/0406455; arXiv;gr-qc/0311019 ; M. Abdalla, S. Nojiri and S. D. Odintsov, Class. Quant. Grav. 22, L35 
(2005) arXiv:hep-th/0409177 ; G. Cognola, E. Elizalde, S. Nojiri, S D. Odintsov and S. Zerbini, JCAP 0502, 010 (2005) 
pirXiv:hep-th/0501096 ; S. Capozziello, V. Cardone and A. Troisi, [arXiv:astro-ph/0501426 G. AUemandi, A. Borowiec, 
M. Francaviglia and S. D. Odintsov, arXiv:gr-qc/0504057; G. AUemandi, M. Francaviglia, M. Ruggiero and A. Tartaglia, 
arXiv:gr-qc/0506123; T. Muhamaki and I. Vilja, arXiv:astro-ph/0506692 ; J. A. R. Cembranos, Phys. Rev. D 73, 
064029 (2006) arXiv:gr-qc/0507039 ; T. Koivisto and H. Kurki-Suonio, arXiv;astro-ph/0509422 , T. Clifton and J. Bar- 
row, arXiv:gr-qc/0509059 , O. Mena, J. Santiago and J. Weller, arXiv:astro-ph/0510453 , M. Amarzguioui, O. Elgaroy, 
D. Mota and T. Muhamaki, arXiv:astro-ph/0510519 : I. Brevik, arXiv:gr-qc/0601100 ; R. Woodard, arXiv:astro-ph/0601672j 
T. Koivisto, arXiv:astro-ph/0602031; S. Perez Bergliaffa, arXiv:gr-qc/0608072, T. Faulkner, M. Tegmark, E. Bunn and 



Y. Mao, arXiv:astro-ph/0612569 , G. Cognola, M. Castaldi and S. Zerbini, arXiv:gr-qc/0701138 , S. Capozziello and 

R. Garattini, arXiv:gr-qc/0702075 , 
[7] S. Nojiri, S. D. Odintsov, Phys. Rev. D 74, (2006) 086005 | arXiv:hep-th/0608008| ; |hep-th/0610 164 S. Capozziello, S. Nojiri, 

S. D. Odintsov and A. Troisi, Phys. Lett. B 639, 135 (200 6). 
[8] S. Fay, S. Nesseris and L. Perivolaropoulos, ,arXiv:gr-qc/0703006[ S. Fay, R. Tavakol and S. Tsujikawa, 

|arXiv:astro-ph/0701479] 

[9] S. Nojiri, S. D. Odintsov and M. Sasaki, Phys. Rev. D 71, 123509 (2005) 'arXiv:hep-th/0504052l. 

[10] M. Sami, A. Toporensky, P. Tretjakov and S. Tsujikawa, Phys. Lett. B 619, 193 (2005) arXiv:hep-th/0504154) ; G. Calcagni, 
S. Tsujikawa and M. Sami, Class. Quant. Grav. 22, 3977 (2005) arXiv:hep-th/0505193 ; S. Tsujikawa and M. Sami, 
arXiv:hep-th/0608178; Z. Guo, N. Ohta and S. Tsujikawa, arXiv:h ep-th/0610336 A. Sanyal, arXiv:astro-ph/0608104 
[11] B. Leith and I. Neupane, arXiv:hep-th/0702002 , B. Carter a nd L Neu pane, JCAP 0606, 004 (2006) arXiv:hep-th/0512262|. 



[12] T. Koivisto and D. Mota, ^arXiv:hep-th/0609155t |arXiv:astro-ph/0606078, T. Sotiriou and E. Barausse, 

|arXiv:gr-qc/0612065 

[13] STTJojiri, S. D. Odintsov and M. Sami, Phys. Rev. D 74, 046004 (2006) 'arXiv:hep-th/0605039 ; D. Konikowska and 
M. Olechowski, arXiv:hep-th/0704.1234; E. Elizalde, S. Jhingan, S. Nojiri, S.D. Odintsov, M. Sami and I. Thongkool, 
arXiv;0705.1211'hep-th]. 

[14] S. Nojiri and S. D. Odintsov, ,arXiv:hep-th/0611071| 

[15] G. Cognola, E. Elizalde, S. Nojhi, S. D. Odintsov and S. Zerbini, [ajXiv:hep-t h/06fll98 

[16] I. Antoniadis, J. Rizos and K. Tamvakis, Nucl. Phys. B 415, 497 (1994); P. Kanti, J. Rizos and K. Tamvakis, Phys. Rev. 

D 59, 083512 (1999); N. Mavromatos and J. Rizos, Phys. Rev. D 62, 124004 (2000); S. Kawai, M. Sakagami and J. Soda, 

Phys. Lett. B 437, 284 (1998); S. Kawai and J. Soda, Phys. Rev. D 59, 063506 (1999). 
[17] S. Nojiri and S. D. Odintsov, Phys. Lett. B 631, 1 (2005) arXiv:hep-th/0508049 ; G. Cognola, E. Elizalde, S. Nojiri, 

S. D. Odintsov and S. Zerbini, Phys. Rev. D 73, 084007 (2006) [a:rXiv:he p-th/06 01C)08| . 
[18] D. N. Spergel et al, arXiv:astro-ph.0603449. 



9 



[19] S. Nojiri and S. D. Odintsov, Phys. Lett. B 599, 137 (2004) ' arXiv:astro-ph /040 3622' ; T. Inagaki, S. Nojiri and 

S. D. Odintsov, JCAP 0506, 010 (2005) arXiv:gr-qc/ 0504054 . 
[20] A. L. Erickcek, T. L. Smit h, M. Kamionkowski, Phys. Rev. D 74, 121501 (2006) [arXiv: astroph /06 10483] ; T. Chiba, 

T. L. Smith, A. L. Erickcek, 'arXiv:astr o-ph/0611867| 
[21] W. Hu and I. Sawicki. arXiv:0 795.1158 



